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1. Introduction 
Let us consider the Hartree equation in Mf^: 

iedt'^'ix, t) = -^A*-'(x, t) + {V{x, t) + U{x, t)) ^^(x, t), 

where 

V{x,t)= [ <Pi\x-y\Wiy,t)\'dy (2) 



1) 



is a self-consistent potential given by a smooth two-body interaction, : M M, even, and 
[/(•, t) : M*^ R for alH > 0, is a smooth external potential (see the next section for the precise 
assumptions on (f) and U). 



In a recent paper [T] the authors of the present one considered the semiclassical limit of the 
version of the Hartree equation corresponding to mixed states, for initial data whose Wigner 
functions do not concentrate at the classical limit. 
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2 SEMICLASSICAL PROPAGATION OF COHERENT STATES FOR THE HARTREE EQUATION 

The problem we deal with in the present paper is the semiclassical asymptotics for ([T]) 
when the initial state is a coherent state centered around the point q,p of the classical phase 
space, namely: 

nix) = e-ia, := (3) 

This problem was studied in [9] in the kinetic (Wigner) picture, see Theoreme IV. 2 therein. 
There it is shown that, under appropriate conditions, the solution IV^ of the Wigner equation 
corresponding to the dynamics ([I]) namely 

i 

e{27ry 
i 



d,W^ + k-d.W^ = j e'^y(v{x + '-^,t)-V{x-'-^,t))dyW%x,k-Od^ + 

J j e'^y{u{x+^-^,t)-U{x~^j-,t))dyW'{x,k-i)di, (4) 



where V{x,t) is the same as in ([T]) equivalently written as 



V{x, t)= U{\x- y\)W'{y, k, t)dkdy, (5) 



(6) 



converges, in weak*-sense, to the solution of the (classical) Vlasov equation 
dtf + k- dj - d^Voix, t) ■ dkf - d,U{x, t) ■ dkf = 0, 



f{x,k,t)\t=o = fo{x,k), 

where 

Vo{x,t) = J (j){\x - y\)f{y,k,t)dkdy, 
and U{x,t) is the same as in ([T]) . The initial condition for (jG]) is given hj Jq = w —* limP^Q. 

e— s>0 

It is easy to check that the conditions of Theoreme IV. 2 in [9] are satisfied for W^{x,v) = 
W'^[^q]{x,v), \E'q as in equation ([3]). In that case (under appropriate assumptions on the pair- 
interaction potential (p and the external potential U) it can be seen that the Wigner measure 
of the wave function verifies 

W'[^']{x,k,t) 6{x-X{t))6{k-K{t)), as £ ^ 0, 

where 

X{t) = Kit), kit) = -VUiXit),t), X(0) = q, KiO) = p. 

In that sense, the semiclassical limit of the problem ([T]) is known to be the Vlasov dynamics 
since it is easy to recognize that, due to the smoothness of the potentials, the limiting measure 
S{x — X{t))6{k — K{t)) is the unique (weak) solution of the Vlasov equation with initial datum 
S{x — q)Sik — p). 

The goal of the present work is to strengthen this approximation. First of all, we construct 
approximations, as opposed to the with weak*-limit, and this yields an explicit control of the 
error in e which allows to recover the shape with which concentrates to a 5 in phase-space. 



semiclassical propagation of coherent states for the hartree equation 

2. Main result 
We will consider the Hartree equation in W^: 

iedt'^^ix, t) = -f A^=(x, t) + {V{x, t) + U{x, t)) ^^(x, t), 

^"(x,0) = %{x), 

where 

V{x,t) = j 0(|x-y|)|vl/-(2/,t)|2rfy (8) 
The initial condition will be of the form 



ao 
qp 



Assumption 2. 

CfM) 3 </> even 



(7) 



and we will make the following assumptions on Oq, and U: 
Assumption 1. 

II«o|Il2 = ||^"oIIl2 = 1, 

x^d^aoix) e for any pair A, B e N'^ with \A\ + \B\ ^3, 

Xi\ao{x)\'^dx = 0, \/i = l...d. (9) 
ki\aQ{k)\^dk = Q, \/i = l...d. (10) 



Assumption 3. 

UeC' (M+,C,3(Mf)). 

Here and henceforth we denote by C^{W^) the space of continuous and uniformly bounded 
functions on M"^ whose all derivatives up the order k are also continuous and uniformly bounded. 



Theorem 2.1. Under Assumptions\^ [H and\^ there exists a constant C such that, V t > 0, 

\m:t) - e^'^^-(*V,V,)IU. ^ C^e^"'^" ■ ^e. (11) 
where Pt is the solution of 

^ (-1 . ^ . <--^'y-"-> ) (12) 

(3o{x) = ao(x), (13) 
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lit) = fj V'mr^)?dr,ds, (14) 

2 

{q{t),p(t)) is the Hamiltonian flow associated with ^ + ^(^'5^) + </'(0) issued from {q,p), 

Cit) := f {p{sf/2 - U{q{s), s) - 0(0)) ds 
Jo 

(the Lagrangian action along such Hamiltonian flow) . 

Remarks: 

• As shown in the proof of the Theorem, the constant C depends only on d, \\U\\w3,°°, 
||0||iy3,oo and sup | |x"^(9^ao| |i,2. 

\A\ + \B\^3 

• Note that in the classical flow the nonlinear potential enters only via the inessential 
constant 0(0). Indeed, due to the symmetry and smoothness of 0, we have 0'(O) = so 
that, in case of concentration as e — )■ 0, the self-consistent field W vanishes. 

• A similar problem for 0"(O) > has been faced in [2] in a semirigorous way. Here 
we treat the case 0"(O) < as well and present an explicit control of momenta and 
derivatives of the solutions (see Lemma 2.3 below) which allow us to estimate the error 
in L\ 

• For a related result (Gross-Pitaevskii equation with a different scaling) see [3]. 

• Assumption [T] can be relaxed by dismissing equation . Indeed even if does not 
hold one can always make a change of variables x ^ x — J x\ao{x)\'^dx. However in that 
case one would have to adjust appropriately the external potential, which of course is 
not translation invariant. 

3. Proofs 

3.1. A Lemma. We first prove the following 

Lemma 3.1. bt{x) := e*'''*^*^ /3i (x) as defined by ![Wil3'f^l4^ is the unique solution of the equation: 
{td, + iA) h{x) = m / |x - v\'\btiv)\'dv hix) + <£:^!M^6,(a;), 

(15) 

bo{x) = ao{x). 

Proof. 

idtbt{x) = -i'{t)bt{x) + e'^^'^idMx). (16) 
By virtue of equations f|T2|) . f|T3l) and f|T^ we find 

idA{x) = ^ / r/2|/3,(r^)|2rf^ 6i(x) + e^^W (-f A(x) + ^x' A(x) + <-.v^c/(.(t),t).> ^^^^^ j ^ 



bQ{x) = ao{x), 



(17) 
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namely 

^^A{x) = -f6,(x) + fflx^ bt{x) + m/ry2|^^(^)|2^^ + <.yMm,t).> ^^^^^^ 

(18) 

bo{x) = ao{x). 

We first notice that the equation (fT2l) for (3t{x) is a linear Schrodinger equation with an 
harmonic potential; therefore the solution /3t(x) of the initial value problem (fT2i) - (fT3l) is uniquely 
determined in L'^{M.'^) and 

IIA|Il2 = ||aolL2 = 1, Vt gM. (19) 
As a consequence of that, it turns out that equation ( |T8l) can be rewritten as 

idAix) = -^bt{x) + m/x2|A(r/)pd7^ btix) + m J^2|^^(^)|2^^ + <.,V^^(gW..).> ^^^^^ 

6o(x) = ao(x). 

(20) 

Furthermore, it is easy to check that if 

xaoix), d.aoix) G L^R''), (21) 

then 

x(3t{x), d,(3t{x) G L2(M^), for all t, (22) 

(see Observation 14.31 below) . 

Condition fl2T|) is satisfied under Assumption [1], so the property fl22|) holds and, in particular, 
there exists a constant C finite for any time t such that 

y \7]\^\Ptiv)\''dv<C, V tGM. (23) 

Thus, by virtue of fl23l) and of Assumptions [H [2] and [3], it follows that the initial value problem 
(|20|) is guaranteed to have a unique solution in and, clearly, ||fcf||j;^2 = Haollia = 1, V t. In 
fact, the equation for bt{x) has turned to be a linear Schrodinger equation with an harmonic 
potential (and all constants appearing in the potential terms are finite thanks to Assumptions 
Eland [3] and to ([23])). 

Now, it remains only to recognize that (!20|) is exactly the same as (fT5|) . To this end it 
is sufficient to observe that, since the equation (fT2|) for /3i(x) is a linear Schrodinger equation 
with an harmonic potential and conditions (Q and (ITOl) are satisfied at time t = 0, we are 
guaranteed that 

%|A(r/)prfr/ = 0, V t. (24) 



Thus, by virtue of (124]) . it follows straightforwardly that (120!) can be rewritten as 

^^Aix) = -f 6,(x) + m / |x - vmm'dv h{x) + <-.v^^(a(t),t).> ^^^^^^ 

(25) 

6o(a;) = ao(x). 
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Finally, it is clear, by the definition of bt{x), that |/3t(x)| = for any x and t. Therefore 

turns to be exactly the same as ffT5]l . □ 

3.2. Proof of Theorem 12.11 We seek an approximate solution to equation ([1]) of the form as 
e.g. in [6l[7l[8l[IIl[T2] 



re/ ,X-q{t) Mt).(.-q(t)) C(t) 

-^''{x,t)=e ■4a( p^,t)e' ^ e' ^ 



where 



q{t)=v{t), p{t) = -VU{q{t),t). 
By inserting the ansatz fl2^ in equation ([T]) we get 



(26) 
(27) 



iedt^^x.t) = e * 



^e^,a{^^^, t) - t^Va{^^-^, t) ■ q{ty 



. p{t)-(x-q{t)) .C{t) 

xe = . 



and 



- C'{t)a{^^,t) - {p{t){x - q{t)) -p{t)q{t))a{^-^,t) 



£^ A ,x — q(t) , p'^(t) ,x — qit) , 



(28) 



_i.A^^(x,t) = e'i 



-i^/eVa{^—^,t) ■ p{t) 



. p(t)-{x-q(t)) C{t) 

e ^ e <^ . 



(29) 



while, with regard to the potential terms in ([T]), we find 

{V{x,t)^U{x,t))W{x,t) = ^-'/'(/ 4>{\x-y\)e-''/'\a{y^^,t)\'dy + U{x,t) 

x-q{t) . p(f).(^-9(t)) .£(t) 
X a( -p-^.Vjd ^ e <^ . (30) 

By (!28|) . (!29|) and (!30|) we get that the amplitude a solves the following initial value problem: 

{%dt + |A) a(/i, t) = i\4(/i, t)a(/x, t) + 

+ i [[/(g(t) + y^/i, t) - Uiq{t),t) - ^VUiq{t),t) ■ fi] a(/i, t), (31) 

a(/i, 0) = ao(/u). 



where 



K(/i,t)= / (0(v^|^-r/|)-0(O))|a(r/,t)|2dr/, 



(32) 
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q{t),p{t) are as in the claim of Theorem 12.11 and we have used the rescahng /i = . 
Note that we should have 

Veifx,t)= [ 0(v^|/i-r/|) |a(r/,t)|2rfr/-0(O), (33) 



instead of (132!) in equation (13T]) . However equation (!3T|) with potential ( !33|l is an Hartree 
equation which preserves the norm so that we can replace (l33l) by (132!) . 
Since G C|(M) is even 0'(O) = and the Taylor expansion yields 



0(v^l/i - ^1) - m = ^0"(O) + eiR{\f^ - r/l), 

(34) 

|i?(|^-r/|)| ^C||niLoo|/i-r/|3, 
while for the terms involving U we find: 

Uiqit) + v^/i, t) - f/(g(t), t) - v^Vt/(g(t), t) ■ /i = e < /i, ^^^M^/x > +£ii?f;(/i, t), 

(35) 

|i?c/(/i,t)KCsup„eN^,l„l=3|V"t/(g(t),t)||/i|3, 

where := V ® V. 

3 3 

The core of the proof is to estimates the two remainders e2R(^\fi — ri\) and e'^Ru{^,t) so 
that we can substitute (0(v^|/i - r/l) - </'(0)) by ^^(/."(O) (as in flMj) ) and f/ (g(t) + v^/i, t) — 
U{q{t),t) - ^VU{q{t),t) ■ /i by £ < /i, ^'^^'^'^V > (as in ([35])). 

In the framework of semiclassical approximation for the linear Schrodinger equation using 
coherent states the method is standard (see e.g. [6], [71 [H [TTl [12]), however we establish these 
estimates again for completeness. 

Denote a((/i) := a{fi,t) and 

htifi) = btiiJ.) - atin). (36) 
By straightforward substitution we get that ho{fi) = (see f[T5|l ) and 

/ \ 



1 ^ /0"(O) r ,2,, , ,,2 , V^U{q{t),t) 



2 



'^"^^^ ^|/i-r/r(|6,(r/)r-|a,(r^)P)rfr/a,(/i) + 



|2, 



i?(|/i - ?7|)|at(?7)| dr] at{fi) -y/eRu{fi,t)at{fi). (37) 

v ' 
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By standard manipulations it turns out that 



d 



Moreover, the term involving Ii can be estimated as follows 



I (Ji, ht) I + (/2, ht) I + v^l {Ru{; t)au ht) \ . 



(38) 



\{h,h)\ < 





< 



< 



l/i - r]\^\ht{r])\ {\bt{r])\ + \at{r])\) dr] at{fi)ht{fi)dfi 



/I ri 



< 



< 2 



while, thanks to flMjl . the term involving I2 is estimated by: 



\{i2,h)\<c\\r 



/ f \fi - r]?\(^t{v)?dr] at{fi)h{fi,t)d^ 



< 



+3(/rfr/|r/P|a,(r/)P)^/^||/z,||^.+3(/d/x|/xha,(/x)p)^/' {J dv\v\\ativ)\') \\ht\\L^ + 

+ (/ci/x|/ina,(/.)p)^/^||a,||i,||/.,|U.). 
One should observe here that / dr] \r]\ \at{r])\'^ < (J dr] (1 + |?7p) |a((?7)p). 
Finally, due to (!35|) . the term involving Ru{fi,t) is controlled as follows: 
\{Ru{;t)a„ht)\ < C sup |V"f/(g(t),t)| f /rf/il/f |ai(/i)||/ii(/i)A < 

a:|a|=3 \J J 

1/2 



< Csup sup |V"f/(g(t),t)| ( / rf/i|/i|^|at(/i)|' 

t a:|Q:|=3 \J 



t\\L^- 



(39) 



(40) 



(41) 



Making use of Lemma 14.21 and equation (!70l) below to estimate terms of the form 



^at\\L2 = (/ rfr/|r/|2™|at(r/)|2)^/^, for m < 3, and 



U dr]\r]\'^\b,ir])\y , for 
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m < 2, in terms of the same quantities evaluated at time t = 0, we easily show, by summing up 
the previous estimates, that there exist three £- independent functions Ci (t) , C2 (t) such that: 

Jt\ML^< V^C, (t ) + C2 (t) 1 1 /li I . (42) 

In particular Ci(t),C2(t) depend on the potentials and U and on the L^-norm of moments 
and derivatives of Qq (up to the order 3). With regard to the time dependence, Ci(t), C2(t) are 
double exponentials Ce'"'^^\ following Lemma [4.21 and observations 14.3^ 14.41 
The conclusion follows with application of the Gronwall lemma. 

□ 

4. Auxiliary results 

Observation 4.1. Observe that under our assumptions the nonlinear equation f31\} can be 
shown to have, for any T > 0, a unique solution in ([0, T], /.^(M'^)) (see e.g. Therefore 
it follows (see e.g. [I3]j that the corresponding time-dependent linear problem 

{idt + |A) u{fi,t) = i / (0(v^|^ - r/l) - 0(0)) \a{r],t)\^dri n(/i,t) + 

+ i {U{q{t) + V^/i, t) - U{q{t),t) - ^VU{q{t),t) ■ ft) u{fi, t), (43) 

m(x, 0) = Mo(x), uoeL'^{M'^) \\uo\\^2 = l, 
has a unique and well-defined propagator. 

Lemma 4.2 (Propagation of Moments and derivatives for a{x,t)). Let a{x,t) be the solution 
of the initial value problem (3l\) . Suppose that for some m G N there exists an e-independent 
constant Mm > such that 

\\x^d^ao\\L2^Mm (44) 

for all A,B eN'^ such that \A\ + \B\ ^ m. 

Form > 2, assume G CJ^{R'^) and U G Ci(Mi+, C™(Mf )). Then, there exists a (finite) 
e-independent constant Cm such that 

llx^Sf a(t)|U. ^ C„e^-^^'"'M^, (45) 

for all A, i? G N'^ such that \A\ + \B\ ^ m. 

Form = 1 inequality holds by assuming G CliW^) and U G C^(]R^, C|(M^)), while 
in the case m = formula ( fT^p becomes an equality and holds with unitary constant (for all t) 
by simply assuming G cf{M.'^) and U G C\Rt , O^t)) ■ 

Remark:The proof makes no use of an energy conservation argument, and this is the reason 
why the Lemma can be established for both signs of 0"(O). 



Proof. Denote 



^^''^{x,t) =x^d^a{x,t), 



(46) 
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e.g. ip^''^{x,t) :=a(x, t). 

It is straightforward to check that 

tdt + I A - '-V,{x, t) - W{q{t) + ^x, t) + \U{q{t),t) + -LV[/(g(t), t) ■ x) t) = 

k=l ^ 



0^«<Afc=l 



+1 E n (t)^r'f/(9W + v^^,^)v^''^(^,t)-^ E n (t)5i(vf/(g(t),t)-x)^^-''^(x,t) 

0</<Afc=l 0<«^A, fc=l 

|«|=1 

(47) 

where B = {Bi, B^, B^) J = {h, ... Jk, ■■■ Jd), A = {A^ Ak, A^) andO < I < A 
means that < < A^ for any k = 1,2, ... ,d. The consistent initial data for fHTj) are defined 

by 

and in particular ip^'^{x,0) := aQ^x). 

Some remarks with regard to our notation are in order; it is clear for example that if 
Bk = or Bk = 1, then the first term on the right-hand side yields no contribution. Similarly 
for Sfc = in the second term and |y4| = for the remaining terms respectively. 
The derivation of P7|) is straightforward by induction. 

Denote by P(t, r) the propagator associated with the left-hand side of equation P7|) . which 
is known to be uniquely well defined in LP' (see Observation 14 . 1 p . As a consequence, for m = 0, 
the result claimed by Lemma [4. 21 follows from the existence of the propagator. We will proceed 
for m e N by induction. 

We will work with vectors including all the moments and derivatives, namely, = 

{i'^"^} A,B:\A\+\B\<m. ^ and 



Os;|A| + |B|s:m 



L2, 



where Xq := L'^{W^). 
For m = 1 we have 



(idt + - -Ve{x,t) - -U{q{t) + ^x,t) + -U{q{t),t) + ^VU{q{t),t) ■ x] 
\^ Z E E E \ ^ J 



-a,^K(a:,t)^°'°(x,t) + -d,U{q{t) + ^ex,t)i,'^'\x,t) - ^d,V{z,t)l=^^,^ ^P''%x,t), 



(4J 
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and 



(idt + - -V,ix,t) - -Uiqit) + ^x,t) + -Uiqit),t) + ^Vf/(g(t),t) ■ x] = 

\ Z E E £ \l E J 

= 7/;"^'°(x,t) (49) 



By virtue of the Duhamel formula we get: 



,\,^=^'\x,t) = P{t,0)ij^^^"{x,0)+ / drP{t,T) 



■dMx,T)ij'''%X,T) 



+ / drP{t,T) 



E ' 



and 

V'°''^(x,t) = P(t,0)7/'°'^^(a;,0) + / dTP{t,r) [4j^^'%x,t 

Jo 

Then, by recalhng that P(t, r) is L^-norm preserving, we find 



(50) 



(51) 



r^'\t)l, < |U--0(0)|U+ / dr 



lL2 



L2 



-d,y,{x,T)ij''^{T) 

E 



+ 



L2 



+ £dr Q9,,f/(g(r) + v/ix,r)--^9.^,t/(g(r),r)^ 



0,0/ 



L2 



and 



||^°''^^(t)L. < ||^°''^^(o)||^,+ / dT\\r^'\T)\\^, 



Now, taking into account the terms involving U in fl52|) . we get 

-5x,t/(g(r) + v^x,r) - ^ 9,^.f/(g(r), r) = 



(52) 



(53) 



diU{z,T) 



\z=q[t) + VSci 



Xj, for some S G (0,e), 



(54) 
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therefore 



L2 



L2 

\ z=q{T)+\/~E X 



< 



L2 



< sup \\d'U{-,T)\\^^\\tP^''^{r)\ 

T£[0,t] 



L2 



(55) 



On the other side, with regard to the term involving Ve in fl52l) . we have 



-5x,V;(x,r) 



< 



< dr] 



</)'(v^|a;-r/|) 



|^°'°(r/,r)r<L / rfr/|x-r/||^"'"(r/,r)r, (56) 



0,0/ 



where L is the global Lipschitz constant of 0' (i.e., the L°°-norm of 0") that is known to be 
finite since G C^iR'^). Then, by ([56]) we get that 



L2 



< dx cir/|x-r/||^°'°(77,r)p / dr/'|x - r/'||^/;°'°(r/', r)n?/^°'°(a;, r)| 



(57) 



< j \x\^\^'''\x,T)\^dx + 3L^ (^J \r]\\^'''\7],T)\^d7]^ ^ 



< L^W W''{r)\\h+3L^\\ krV)lli2, 

where we made use of the fact that ||'0°'°(r)||j;^2 = ||'0°'°(O)||i2 = ||ao||L2 = 1, for any time r. 

At this point we observe that || |x|^°'°(r)|||2 = || |x|a(r)||^2 = II'V^°'^^(''")IIl2- So that, we 

j 

have just proven that there exists a constant C > depending only on the L°°-norm of the 
second derivative of 6, such that 



-9..V;(x,r)^°'°(r 



(58) 



Now, by using fISS]) and f lSS]) in f lS^ . we obtain that 



T 



Il2 ) 



(59) 



where C is not the same constant of formula fl58l) - we denoted it by the same symbol just for 
the sake of simplicity - since here it is depending on 0, as previously, but even on U (through 
the L°°-norm of its second derivative, according to fl55|l ). 

Now after f l55|) . summing over j = 1, . . . , ci in equations ( 159|) and fl53|) and then adding them. 
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we get 

0mx.<0m\x.+CjdT0{T)\\x,. (60) 



The conclusion follows by applying the Gronwall lemma, i.e. 

0mx. < ll^(0)|ke^*<Mie^*, (61) 



where Mi has been defined in 

For m > 2, the previous inductive step from m = 1 applies almost verbatim: first, by 
virtue of the Duhamel formula, we write the solution of equation f H7|) by using the propagator 
P(t, r) associated with the time evolution on the left-hand side. Then, by using the L^-control 
on P(t, r), it only remains to show that the "source terms" appearing on the right-hand side 

of (147|) are suitably uniformly bounded in terms of 1 1?/;^'-^! 1^2 or H^E^Hxm- The way to do that 
is by using 1 1?/^^'^! 1^2, |A| + < m as constants now. 

For example, let us look at the term involving the potential on the right-hand side of 
iZD, i.e. 

\A-l\=l 



^\ E i\[^^^Qt'ys{^M^''{^.t). (62) 

n^;^ A I — 1 \ k / 



0!il<A k=l 
\A-1\>1 

The estimation for any of the terms in the last sum reads as: 

\\dt' I i0( v^k - V\M''%V, t)\'drj V^''^(x, t) I \l, ^ 

^ (£^^||</.(^-')(x)||ioo)'||/|x-r/||V'°'°(r/,t)prfr/ ^''^(x,t)||2, ^ 

<:2D\\ |7^i<'0(t)iu2||^''^(t)iu2|| \x\^p^'^mL^+D\\ \v\^''%mL.\\ij'^^mh+D\\ ixi^p^'^^mh 

(63) 

where D is a constant only depending on H^*-^ ')(x)||2,oo (that is finite under our assumptions 
since A — I < m). Furthermore, it is clear that, by construction, we are guaranteed that the 
exponent for e is non negative. 

On the other side, the estimate for any of the terms in the first sum on the right-hand side of 
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( l62l) is given by 

\\dt' I -MV^\x - vDH'^'^v, t)\'dv t) I li. < 

^ L\\J\x-7]\\iP^'\r],t)\^d7] ip^'^{x,t)\\l, < 
^2L|| |r/|^0'°(t)|U.|l^''^WIlL^II + \v\i;''%t)\\L.\\^^'^{t)\\l, + L\\ \x\^''^ 

(64) 

where L is the global Lipshitz constant of 0' (see (156|) ). which is guaranteed to be finite since 
4) e Cfe™(R'^), with m > 2. 

Now, by virtue of the estimate we proved for m = 1 (see (!60|) ). from (!63|) and we find 

that 



(65) 

L2, 



where K = max{D, L} and and we recall that |/| + < — 1 + < m — 1 and |/| + + 1 < 
\A\ - 1 + \B\ + 1 < m. Thus: 

\\dt' I i0(v^k - V\W'%V, ^''^(x, t) I \l, < K{M,e^' + 1)0 (t) \ (66) 

Concerning the terms involving the potential U on the right-hand side of (1471) . the idea is 
quite similar. In fact, we observe that 

-s E U{'^^)^t'U{q{t) + V^x,t)^P^'^{x,t)- 
Oi:i<Ak=l 

-7s E u{l:)dim{q{t),t).x)r-''^{x,t) = 

\l\=l 

= E CA,i,BUd,Uiqit) + ^x,t)^^-^'^{x,t)-j=^dAVU{qit),t)-x)^'^^^^^ ^^^^ 
\l\=\ 

+7 E X\{'ll)dt'Um'r^exM^''{^.t). 

Os£/<A, fe=l 
\A-l\>\ 

where we made a discrete change of variable I A — l \\y the first term of the left-hand side. 
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Now, with regard to first term of tlie right-liand side, tlie estimation tliat lias to be used 
is exactly the one we did in flSSl) . thus one finds, VZ : |/| = 1 



e Je 



< 



<sup 9^f/(-,t) L U <sup 



(68) 



(the adjustment for / = is obvious). 
Now, for the last term in (|67l) we have 



-^^t'U{q{t) + V^x,t)^p^'^{x,t) 



< sup I |9(^-^)f/(-, t) I WiJ^'^'it) I 

L2 i 



L2 



< 



< 



p(-)[/(-,t)IU.||l^(t)|| 



sup 



(69) 



where we used that A - I < A < m, |A-/|-2>0 and l + B<A + B<m. 

Similar (simpler, in fact) estimates can be shown for the other terms on the right-hand 
side of (HZD. 

□ 



Observation 4.3. [Propagation of moments and derivatives for /3t(x)] /3((x) was defined in 
equations [W\). / f73]) . Under the assumptions of Lemma \4.^ regularity estimates for f3t{x) 
analogous to Lemma \4-S\ for a{x,t) hold, i.e., for any t > 



\A'\ + \B'\^m 



x^d^aollL^, yA,Be N'^: \A\ + \B\ < m. 



Remarks: 

• The proof is in fact simpler with respect to the one of Lemma I4.2t it can be checked 
easily that, due to the fact that we have to deal with harmonic potentials, the terms 
that arise from the differentiation of the potentials turn to be exactly of the form x/3f (x) 

(~ II \E'/3(t)||xi if we denote by 4'f'^{x,t) the quantities x^d^(3t{x) and we define ^'/^(t) 
consistently), i.e., precisely the kind of objects we want to recover to apply the Gronwall 
Lemma (see the proof of Lemma [4. 2p . 

• As a consequence of Observation 14. 3[ by Assumptions [H [21 [3] we are guaranteed that, 
in particular, there exists a e-independent constant C > depending on the L°°-norm 
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of the second x-derivative of U{x,t) and on |0"(O)|, such that 

j dx < (^J dx |x|>o(a;)r^ e^* < oo, VI 

We remind that this is exactly what we need to make the proof of Theorem 12.11 work 
succesfuUy (see (!23|) ). 

Observation 4.4. [Propagation of Moments and derivatives for bt{x)] Although apparently 
bt{x) solves a nonlinear equation, it can be obtained as the solution of a linear Schrddinger 
equation with an harmonic potential whose coefficients are determined by the Li^-norm of the 
first moment of f3t{x), by 0"(O) and V'^U (see / f^) and Lemma l3J\) . 

Therefore, as a consequence of Observation \4^ it follows that, as long asU E (M^^, C™(]Rf )) 
and (j)"{0) is finite, we can get a result for bt{x) e.g. analogous to Lemma \4.S\ for a{x,t), i.e. 

x^d%eL\ yA,Be N^: \A\ + \B\<m V t, 

under the same assumption ( [^^P on the (common) initial datum ao(x). 

Note that, in particular, by Assumptions [H |2l [3] we are guaranteed that there exists a e- 
independent constant Ct > depending on the L°°-norm of the second x-derivative of U{x,t), 
on 10" (0)1 and on time t (but finite for any t), such that 

J dx |x|2™|6t(x)|2 <Ct, V t, m < 3. (70) 

We observe that (!70|) . for m = 2, is exactly what we need to make the proof of Theorem 12.11 
work succesfuUy (see fl39l). fliOl) and fl4Tl)). 



5. Higher order approximations 

On the basis of the above results, it seems natural to ask whether it is possible to go beyond 
the -y/^-approximation discussed previously (see (fTT!) ) and to find higher order corrections a^^\ij) 
to the amplitude af\ij) := bt{jj) so that the right-hand side of f lTT]) gets of size of any power of 
e, as this is the case for the linear Schrodinger equation [TTl [12]. Although we will not present 
all the (tedious) details of the construction, we claim that one can determine a semiclassical 
expansion 

af(/i) = af\ii) + v^af ^(/i) +eaf\ii) + ■■■ + e'/'af\ii) + . . . , (71) 

with 

4''^ (/^) = 4,0^0 (At), (72) 

such that 

vl/^(a;,t) = e^^+^^(*V^%^(^) + 0(6-). 
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In order to determine the equations governing the evolution for each coefficient ai"* (/i) we need 
to look at the expansion for the potential terms appearing in (l3Til . With regard to the nonlinear 
part involving the pair interaction 0, we get: 

\{^{V~e\^i-v\)-m) = ^-^<P'iO) + ^^r(0) + ^'^3~^'' 0"(O) + 

+ ■■■+ (v^)'l/;-^l% W(o) + ... (73) 

In Theorem 12.11 we were assuming G C^(M). Clearly, if we want to go to higher orders in 
the approximation we need more smoothness on (p and on the external potential U. Therefore, 
here and henceforth we assume: 

0gC6~(M), and even (74) 

so that we have 

] (0(v^l/i - ^1) - 0(0)) = ^^^0"(O) + ■ ■ ■ + (v^)^"-^^^^^0(^"Ho) + . . . n > 2 (75) 

Observation 5.1. Assumption f/4\ ) is actually too strong if one wants to deal with an approx- 
imation up to a certain order k. 

With regard to the linear terms in ( 1311) involving the external potential U, we get: 

i m<lit) + y^/i, t) - U{q{t),t) - ^VU{q{t),t) ■ fi) = 

=< /i, ^^^/i >+■■■ + (v/i)-^^:^^ 
where of course n ^ 3. Here we are using the notation 

Analogously to what we observed for the pair-interaction 0, we need more smoothness for 
U, so that here and henceforth we require: 

UeC\Rt,C^iRi)). (77) 

Now, inserting ( I7T]) . (175|) and ( !76|) in (13T]) we readily arrive to a sequence of problems for 
the coefficients a^^\^) of the expansion fl7T]) . For /c = we obviously find: 

4m, 0^(0) /',^,„ ^,2,J0),„„2, V2t/(g(t),t) 
2 2 

(0) 



+ + ^ / - ri?\af\r,)\^+ < " "T'' V > ) 4° V) = 0. ^^^^ 
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namely, the initial value problem that we had for bt{fi) in the previous sections (see ([I5 
for k = 1, we find: 

V'U{q{t),t) 



Then, 



0"(O) 



4'^(/i) = 0. 



r(o) 



dri\fi-r]\'^\a^^\ri)\'^+ < fi 



-/i > a 



.(1) 



dr^\f^-r^\'2^af^\^)a\'>iv)]]cir{fi) + 



V'^(9(t),t) 
3! 



(79) 



This is a linear initial value problem where the left-hand side is known to have a unique well- 
defined L^-propagator P^'^^ (t) due to the existence and uniqueness in of the solution af'^ (/i) 
of the zero-order initial value problem ( |78l) and to the L^-control on its first moment (see 
Observation 14.41 and flTOj) ). Then, it is easy to see that, writing the solution a[^\iJ,) through 
the Duhamel formula (with "leading" propagator P^^\t)), the well-posedeness in for fl79l) is 



jj" af^ (/i) (which is achieved thanks 



guaranteed by the L^-control on the source term 

to the smoothness of U and to the L^-control on the third moment of af\n) - see (ITOil ) and 
by the following estimate: 

y'(o) 



,(0) 



L2 



+ 



2 



d^ 



dfi 



drilfi — r] 
drilfi 



+ 



< C 



(0) 



,(1) 



L2 



(80) 



Here C^(o) > is a constant that only depends on the moments of af'\jj,) up to order 3, 
depending on the potentials U (trought the L°°-norm of its second derivative) and (through 
the quantity |0"(O)|) and on the initial derivatives and moments of af'\fi) up to the order 3 (see 
Observation 14.41) . that, as in the previous sections, we assume to be finite. By virtue of (l80!l 
and the L^-control on the term involving U in (!79|) . the Duhamel formula and the Gronwall 
lemma allow to conclude that 



,(1) 



V t 



Moreover, following the same lines of the proofs presented and discussed in the previous sections 
(see Lemma 14.21 Observation 14.41 and subsequent remarks), it can be easily checked that by 
assuming enough regularity for the (zero-order) "full" initial datum a'"Q\^) = ao(/i), in such 
a way that we control in a sufficiently high number of moments and derivatives, we can 
control the derivatives and moments of a[^\jj,) up to any fixed order m, i.e: 

x'^d^ai^^ G L2(M^), V t, V A, 5 e N*^ : |A| + |5| < m. (82) 
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This will be crucial to go on with the higher orders dynamics because, for example, the equation 
for the second coefficient af^ {fj) is 

+ ^ (/ 'ir, b - „p2S|a™(,)al"(„)l) »<"(„) + . .^a^ip). 

(83) 

So, again, as for the case k = 1, we obtained a linear initial value problem where the propagator 
associated with the left-hand side is P^'^\t), that is known to be uniquely well-defined in L^. 
Then, as before, the solution a), (/i) can be written through the Duhamel formula, applying the 
propagator P(°)(t) to the term (^J dri\fi — ri\'^2^\af'\r])a[^\r])]j af\^) and to the various 

source terms in ( !83l) . The term which is linear in a\ (/i) is estimated as in ( !80l) while the source 
terms are controlled in by virtue of the control on moments and derivatives of o!i\^) and 
a|^^(/i). In the end, by using the Gronwall lemma, we get 

af^ G L^{W), yt. (84) 

and, moreover, by assuming a sufficiently high number of moments and derivatives of the (zero- 
order) "full" initial datum ao°''(/x) = ao(/x) to be controlled in L^, we can control as well the 
derivatives and moments of a), (/i) up to any fixed order m, i.e: 

a:^af af ^ G ^^(M^), V t, ^ A, B e : \A\ + \B\ < m. (85) 

At this point it is clear how to proceed in general. The equation for af^^^) is a linear 
Schrodinger equation with a source term involving the coefficients aj"''(/i) with n < k, which 
have been estimated by the previous steps. The L^-control of a^*^'' (/i) follows by the L^-control 
on a sufficiently high number of moments and derivatives of a["'\fi) with n < k. 

Acknowledgments. We would like to thank R. Carles for pointing out that an extra hypoth- 
esis was needed, and a mistake in the remainder which appears in the main theorem. 
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